PREDICTION OF THE IN-GAP STATES ABOVE THE TOP OF THE VALENCE 
BAND IN THE UNDOPED INSULATING CUPRATES DUE TO SPIN-POLARON 
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In the framework of the generalized tight binding method we have calculated the quasiparticle 
band structure and the spectral functions of the undoped cuprates like LaiCuO/±, SriCuOzCl-z 
etc. Due to spin fluctuations the in-gap state appears above the top of the valence band in the 
undoped antiferromagnetic insulator similar to in-gap states induced by hole doping. In the ARPES 
experiments the in-gap states can be detected as weak low energy satellites. 

PACS numbers: 74.25.Jb, 71.18.+y, 79.60.-i 



I. INTRODUCTION 

The key issue to understand the nature of high- 
temperature superconductivity in the cuprates is the evo- 
lution of the electronic structure from an antiferromag- 
netic insulator to a superconductor with hole doping. 
The appearance of the in-gap states above the top of the 
valence band in slightly doped cuprates has been found 
experimentally [1-4]. In the metallic underdoped regime 
the ARPES measurements [5] reveal the concentration 
dependent band structure of Bi-2212. With improving 
the ARPES resolution recently the formation of the new 
quasiparticle states at the transition from insulator to 
metal in La2— x Sr x CuOi has been found [6, 7]. 

The formation of the in-gap states with doping has 
been obtained theoretically in the numerical studies of 
small cluster in the framework of the t — J model, Hub- 
bard model, and the 3-band p — d model [8-11]. The band 
structure calculations of the C1/O2 layer in the framework 
of the multiband p — d model by the generalized tight- 
binding (GTB) method [12] with account for strong elec- 
tron correlations have revealed the unusual in-gap state 
at the top of the valence band with zero spectral weight 
for undoped insulator that acquire the dispersion and 
non-zero spectral weight with hole doping [13]. In all 
models of strongly correlated electrons the hopping of 
hole in the antiferromagnetic background is renormalized 
by spin fluctuations. To clarify the origin of the in-gap 
state we have studied in this paper a spin-polaron effect 
both analytically in the framework of the t — t' — J model 
and numerically by the GTB method similar to [13]. We 
have found that the spin fluctuations provides non-zero 
spectral weight and dispersion of the in-gap state as well 
as hole doping does - even without doping the spin ex- 
citations that are present in the antiferromagnetic state 
due to the quantum spin fluctuations at all temperatures 
including T = results in the non-zero in-gap spectral 
weight and dispersion above the top of the valence band. 



This state can be detected by the ARPES measurement 
as a weak satellite at the low energy shoulder of the 
main peak. Moreover the concentration of spin fluctu- 
ations n s f increases with temperature, and we expect 
the growth of the in-gap spectral weight ~ n s f . 



II. GTB METHOD RESULTS 

A dispersion equation of the GTB method for the 
quasiparticle band structure of the CuOi layer looks like 
[13] 
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Here to is a quasiparticle band indexes given by a pair 
(p, q) of the initial and final multielectron configurations 
E p (n + 1) and E q (n), Q m = E p (n + 1) — E q (n) is a local 
excitation energy. The local excitation \q) — > |p) is de- 
scribed by the Hubbard operator X pq — \p) (q\, a filling 
factor F(m) = {Xpp) + Two magnetic sublattices 

are denoted by indexes A and B, and a is a spin projec- 
tion. The interatomic hopping is Tj$, where the single 
hole basis set A includes 5 orbitals: copper d(x 2 — y 2 ), 
d(3z 2 — r 2 ), in-plane oxygen p(x), p(y) and apical oxygen 
p(z), j\cr(m) is a parameter of a single hole annihilation 
operator in term of the Hubbard operators 



a\a = 7a. 



r (m)X r 



(2) 



'Electronic address: sgo@iph.krasn.ru 



The local multielectron energies and parameters 7x0- (m) 
are obtained after the exact diagonalization of the multi- 
band p — d model Hamiltonian for the unit cell. In 
our case the unit cell is CuOi cluster for La-iCuO^ and 
C11O4CI2 for Sr2CuC>2Cl2- The similar equation has 
been known long ago for the nondegenerate Hubbard 
model as the Hubbard I solution and has been used re- 
cently to study magnetic properties of transition metals 
[14, 15]. 
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The essential for undoped cuprates multielectron con- 
figurations are d 10 p 6 (vacuum state |0) in a hole rep- 
resentation), single- hole configurations d 9 p 6 , d 1Q p 5 , and 
two- hole configurations d s p 5 , d 9 p 5 , d 10 p 4 , d 10 p 5 p 5 . The 
minimal energy in the single-hole sector of the Hilbcrt 
space has the b\ g molecular orbital, and in the two-hole 
sector the 1 A\ g singlet that besides Zhang-Rice singlet 
contains several more local singlets. A staggered mag- 
netic field split b\g levels by spin: 

£Aa = s% - crh, e Ba =£i+crh. (3) 

The top of the valence band is given by the quasi- 
particles with to = 1: X\ = |6i 9 ,f) ( lj 4ig | and X\ = 
\b\g.i) A\ g \, as usually there is spin degeneracy of the 
band in the antifcrromagnetic state. The occupation 
number n p = (X pp ) are calculated self-consistently via 
the chemical potential equation. In the mean-field Hub- 
bard I approximation the solution of this equation for the 
hole-doped cuprates with hole concentration = 1 + x 
is given by 

?ii T = n A i(bi g ) = l-x,mi = °' n 2 = n^Aig) = x. (4) 

For the band to = 1 we get Fj^(V) = 1 while for the 
band to = 2 with X\ = |&i g ,x) C^-isI the filling factor is 
Fa|(2) = x. The quasiparticle spectral weight is propor- 
tional to the filling factor, thus it is the band to = 2 that 
forms the in-gap state. In the limit x — > its spectral 
weight is zero, when x ^ this band acquires both dis- 
persion and nonzero spectral weight. The corresponding 
concentration dependent bands structure has been ob- 
tained for La2- x Sr x CuOi in [13] and the chemical po- 
tential fj.(x) dependence, the Fermi surface evolution with 
doping have been studied in [16]. 

To go beyond the mean-field Hubbard I approximation 
one has to calculate single-loop diagrams for the self- 
energy [17]. In the ferromagnetic or antifcrromagnetic 
state the most important contribution is given by loops 
with spin- wave excitations [18] (a spin-polaron effect). 
According to [18], the main effect of the spin excitations 
is given by the spin-wave renormalization of the multi- 
electron configuration's occupation numbers, so instead 
of one gets 

nii = C 1 _ x )( l ~ n *f)> n U = C 1 - x)n sf ,n 2 = x, (5) 

where n s f is the occupation of the spin-minority level and 
it determines the spin-fluctuation decrease of the sublat- 
tice magnetization 

(S z A ) = (l-x)(l/2-n sf ). (6) 

Concentration of the spin fluctuations is equal to 2n s f. 

Thus the filling factors for the valence band F(l) = 
1 — n s f, and for the in-gap states F{2) = x + n s f. It 
means that the spin-polaron effect results in the non- 
zero spectral weight of the in-gap states even for undoped 
cuprates La,2CuO^ and Sr2Cu02Cli- The quasiparticle 
band structure and the spectral function for the undoped 
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FIG. 1: The quasiparticle band structure (a) and the spectral 
function (b) of the undoped La2CuOi with a spin-fluctuations 
n s f = 0.2 calculated by the GTB method. The Fermi level is 
above all bands shown here 



La2CuOi with n s f = 0.1 are given in the Fig.^ Here the 
lowest band is formed by hole hopping via 2-hole triplet 
3 Bi g state - this aspect was discussed in [13, 19]. The 
next band (to = 1) is the top of the valence band with- 
out spin fluctuations with a maximum at k = (7r/2, 7t/2). 
The upper band (to = 2) formed by the dispersion of 
the in-gap states. Despite of its width each state has a 
low spectral weight as seen in the Fig. \T]p and the total 
number of states in this in-gap band (without doping) is 
equal to n s f. The appearance of such non- Fermi liquid 
states is the direct effect of strong electron correlations. 
The maximal spectral weight of the in-gap state is near 
(tt, 0) point of the Brillouin zone (BZ). At the (tt/2, tt/2) 
point the two bands are degenerate, and we cannot sep- 
arate the contribution of the in-gap band to the spectral 
function Ak(E). 
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III. t — t' — J MODEL TREATMENT 



where 



To clarify the properties of the in-gap band we study 
the spin-polaron effect in the t — t' — J model, which is an 
effective low energy model for the multiband p — d model 
[20], with the Hamiltonian: 
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where S / are spin operators and n t are number of parti- 
cle operator, tf g and Jf g are the hopping and exchange 
integrals correspondingly. In the Hubbard 1 approxima- 
tion it is easy to obtain in the undoped case the follow- 
ing intrasublattice and intersublattice Green functions 

(feA,geB) 



((xy\X'f°)) E = ±J2Gl A (E)e 

k 



ik(f-f') 



ik(g-f') 



Matrix Grin function in momentum space could be writ- 
ten as: 
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Jn n aa, with a = A, B. 



{^0 ~ ^k ) n a 

Here t k and t k (Jq and Jq) are the hoppings (ex- 
changes) in momentum space between different and the 
same sublattices correspondingly. In simple case of next- 
nearest-neighbors approximation we have 

t k = 2t (cos k x + cos k y ) , t A = At' cos k x cos k y , 



J B = 4J, 4 = 4J', 

with primed values corresponding to next-nearest hop- 
pings and exchanges. The occupation factors of one par- 
ticle state with different spin projections are denoted by 
riAo and nsa- In the mean field Hubbard I approxima- 
tion riA-\ = (1 — n s f) and uai = at T — 0. Using the 
same arguments as in the Section [HJ we go beyond the 
Hubbard I approximation by renormalization the occu- 
pation numbers with spin fluctuations. It results in the 
undoped La^CuOtx in 



n A\ = (1 - n sf ),n Ai 



(9) 



The condition D = gives two branches of quasiparti- 
cle spectrum: 



E 



fcT 
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(10) 



(3 k = {4 - J S + J A f (1 - 2n s/ ) 2 +4 (if ) 2 (I - n. f ) n sf . 

If we set concentration of the magnons to zero we im- 
mediately get one dispersionless state and one dispersive 
state with dispersion governed by hoppings between dif- 
ferent sublattices: 



E 



— E\ -/X- J , E k * — E\ -/i + t fc - Jq. 

n s f—0 n s f—0 

(11) 

If the values of inter-sublattice hoppings and exchange 
are small then the difference between two energy levels 
is of order J: 



AM 



E 



fcT 



E; 



4J. 



n s f — 



(12) 



In Fig. [21 the quasiparticle dispersions corresponding 
to equations (|I0|l and {TJ arc shown. Parameters were 
taken from effective low-energy model [20] of multiband 
p - d model and equal to t = -0.587, t'/t = -0.085, 
J/t = 0.392, J' jt = 0.00037. The value n sf is calculated 
in the effective quasi-two-dimensional Heisenberg anti- 
ferromagnetic model, n s f = 0.2 for typical in LaiCuO/± 
ratio 10 -5 of the interplane and intraplane exchange pa- 
rameters [21]. 

The distance between two spectrum branches for non- 
zero concentration of the magnons is less then distance 
()12JI for n s f = by factor proportional to (1 — 2n s f) = 
2(5'). 

The lower quasiparticle branch for n s f = 0.2 (but not 
for n s f = 0!) clearly resembles dispersion obtained in 
self-consistent Born approximation [22] and GTB method 
[13]. This proves that two different approaches to treat 
spin fluctuations lead to the similar results. 

Introducing energy difference AE k = E 



"fcT 



E k ^ we can 



write down spectral functions A ka (E) 
in the form: 



—Im 



SpG k 



A k1 (E) = ui5[E-E+)+viS[E-E 



J k\ 
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( j b 
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2AE k 



4 



(13) 
(14) 



"fc- 



Obviously, for n s t = the value u k = 0, v k = 1 and 
there will be the non-zero spectral function A k -\ (E) = 

s(e-e 



fc1 



corresponding to only one dispersive state 

(tTT|l . In Fig. Eh the spectral functions versus energy 
for different symmetric points in momentum space are 
shown. Comparison of spectral intensities in case of pres- 
ence and absence of n s f (solid and dash-dotted lines) 
indicates that the second satellite peak appears above 
the main peak in (tt/2,w/2) and (tt, 0) points. It is the 
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FIG. 2: Quasiparticle dispersion (a) and spectral function 
peaks (b) in the t — J model for x — and n s f = 0.2 (solid 
lines) and n s f = (dot-dashed lines) 



satellite peak that represents the in-gap state. In the 
(7r/2,7r/2) point the distance between two peaks is pro- 
portional to J (see equation (|12H ) but at (ir,0) point 
the distance is proportional to |J + t'| (or, generally, 
I Jq — ijn) and will also take place even at zero J. The 
last statement emphasizes importance of next-nearest- 
neighbor hoppings t' at low doping. Indeed, if one con- 
sider t — J model with only nearest-neighbor hoppings 
then the doped hole will not even be able to move with- 
out spin fluctuations (n s / = 0, see equation (|ll|l ) and 
at low n s f the dispersion will be governed by (t£ — Jq) 
term in equation (|10|l but not by nearest neighbor hop- 
ping if. 

Now we will discuss the higher-order corrections to pre- 
vious results. Two main effects are expected: i) quasi- 
particle decay and ii) renormalization of the real part of 
the self-energy. The main change with introduction of fi- 
nite quasiparticle lifetime will be broadening of spectral 
peaks. This effect indirectly presented in Fig. [3j where 
the delta function peaks in spectral function are arti- 
ficially broadened by the Lorentzian. To analyze renor- 
malization of the real part of self-energy we will use more 
rigorous approximation - generalized Hartree-Fock ap- 



proximation [23]. In this approximation the equation 
of motion for operator X® a is renormalized by two-sites 
static correlation functions: 

ij t Xf = [( £l - M ) + M fa ] Xf + Tf g ,*X /, (15) 

g 

where Mf a and T/ giCr are the renormalized chemical po- 
tential (exchange integral) and hopping integrals corre- 
spondingly: 

Mfa = £i /ff (XfXf) 

g 

- E j / 9 [( x ? x r) + <*r*r > - <*r*r>] . 

r f „ = t Sg [{XfXf) + {XfX?) + {XfXf) 

+ {x?x™) + {xfxf)] + J fg (xfxf) . 

It is clear that the equation H15[) has the same lin- 
earized form as in Hubbard I approximation but with 
renormalized chemical potential and hopping integrals. 
It means the qualitative results of Hubbard I considera- 
tion will be the same but quantitatively they may change. 
Namely, due to renormalization of the exchange inte- 
gral the distance between the in-gap and main spectral 
peaks will be shorter then expected from l|12|l and the 
concentration dependence of peak positions will appear 
due to both renormalizations. Meanwhile, the underly- 
ing physics of the in-gap state will be unchanged and its 
dispersion will be governed by spin fluctuations. 



IV. CONCLUSION 

It is clear from the spectral function both in the Fig.^ 3 
and Fig. [2) that there is a pseudogap between the in- 
gap band and the valence band, for the undoped cuprate 
both bands are occupied and the chemical potential lies 
above the in-gap band. With doping fi(x) is pinned to the 
in-gap state [16] up to optimal doping. The pseudogap 
is fc-dependent. In the (ir/2,ir/2) point of the BZ the 
value of the gap is AE (n/2,n/2) ~ J (1 — 2n s f), while in 
the (tt, 0) point AE (tt, 0) ~ | J + t'| (1 - 2n sf ). We can 
compare results of the p — d model and the t — J model in 
our solution only in the limit U — ► oo (J — > 0), then we 
get AE(n/2,n/2) -> and AE (n, 0) -> |t'|(l-2n s/ ) 
that corresponds to the Fig. At J ^ there is the 
additional contribution to the pseudogap, and we may 
expect the in-gap satellite both in (tt/2,tt/2) and (tt,0) 
points of the BZ. 

In conclusion, we have shown that the spin-polaron 
effect results in the formation of the in-gap band above 
the top of the valence band even in the undoped cuprates. 
The spectral function of the in-gap states has a form of 
small low energy satellite that can be detected by ARPES 
measurements. The most interesting to the ARPES stud- 
ies are (%, 0) and (ir/2, 7r/2) points of the BZ. For the hole 
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doped cuprates there are two contributions to the in-gap 
spectral weight: the mean field contribution is given by 
doping concentration x and the spin-fluctuation contri- 
bution is given by the magnon concentration 2n s f. The 
latter is temperature dependent resulting in the increas- 
ing satellite intensity with the temperature growth. 
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